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1.  Introduction 


Complex  metabolic  networks  contain  many  essential  enzymes  that  cannot  be 
genetically  altered  or  deleted  without  adverse  impact  on  the  host  cell.  We 
envision  an  alternate  strategy  to  modulate  metabolite  flow  through  these  essential 
enzymes,  which  will  enable  unprecedented  metabolic  control,  leading  to  more 
efficient  production  of  commodity  chemicals  (i.e.,  biofuels),  better  waste 
remediation,  and  novel  antibiotics.  Our  2  goals  are  to  computationally  discover 
and  design  small  molecules  and  modified  bacterial  genes  (for  producing  peptide 
tails  to  a  common  bacterial  protein)  that  will  inhibit  native  bacterial  metabolic 
enzymes  and  alter  metabolic  output.  In  this  one-year  effort,  we  attempted  to 
inhibit,  via  small  molecules,  the  essential  hydrogenase  (HydA)  (Fig.  1)  of 
Clostridium  acetobutylicum  (C.  acet.),  ’  which  is  predicted  to  subsequently 

•5 

increase  butanol  solvent  production  by  the  organism.  Unmodulated  wild-type 
C.  acet.  yields  a  3:6:1  ratio  of  acetone,  butanol,  and  ethanol.  External  modulators 
in  the  form  of  small  molecules  would  provide  a  proof  of  principle  that  we  can 
inhibit  the  HydA  and  shift  metabolism  toward  butanol  production. 

Inhibitory  Active  site 


Fig.  1  Hypothetical  internal  peptidic  inhibition  of  bacterial  HydA 

This  report  provides  a  structural  template  for  a  future  study:  the  design  of  small 
peptides  to  inhibit  in  the  active  site  of  the  HydA.  This  peptide  sequence,  in  turn, 
may  be  inserted  at  the  head  or  tail  of  the  sequence  of  an  abundant  bacterial  protein 
to  result  in  an  engineered  inhibitor  of  the  HydA.  The  benefit  of  this  approach  is 
that  the  concentration  of  the  inhibitory  peptide  will  be  directly  tied  to  the  presence 
of  a  commonly  expressed  bacterial  protein  (CEBP).  As  a  bacterium  undergoes 
stress  from  reduced  HydA  activity,  expression  of  the  CEBP  will  drop,  which  in 
turn  will  reduce  concentration  of  the  designed  inhibitory  peptide,  thus  creating  a 
feedback  loop  that  will  limit  cell  stress  and  death. 
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2.  Objectives 


The  objectives  in  this  report  are  the  following: 

•  Perfonn  a  computational  screen  of  one  million  compounds  against  2 
distinct  surface  regions  of  the  C.  acet.  HydA  protein  structure  (based  on 
Peters  et  al  4)  and  select  a  total  of  20  compounds  for  experimental  testing. 

•  Using  experiment,  determine  if  any  of  the  20  small  molecules  inhibit  the 
bacterial  HydA,  thereby  reducing  hydrogen  output. 

The  following  are  future  objectives: 

•  Analyze  chemical  functionality  and  predicted  binding  mode  of  all  of  the 
small  molecule  hits  to  design  small  peptides  that  will  inhibit  the  active 
site. 

•  Genetically  engineer  the  top  peptide  candidates  into  the  C-terminus 
(or  N-terminus)  of  an  abundant  cytoplasmic  protein. 

•  Test  the  bacterial  system  for  viability  and  hydrogen  production. 

Successful  completion  of  these  goals  would  be  a  strong  scientific  result,  yielding 
future  research  directions  for  the  US  Anny  Research  Laboratory  (ARL) 
biotechnology  program  and  strengthening  the  repertoire  of  ARL’s  Computational 
and  Information  Sciences  Directorate  and  Sensors  and  Electron  Devices 
Directorate  (SEDD)  scientific  tools  for  collaboration  with  life  and  biomaterials 
scientists  within  the  Department  of  Defense  and  external  institutions. 

3.  Methods 


3.1  Computational  Modeling  of  the  Protein  and  Ligand 

We  used  various  modeling  tools  to  determine  the  correct  site  for  screening 
inhibitors  of  the  enzyme.  The  methods  we  used  are  described  in  detail  in  Section 
4  of  this  report. 

3.2  Computational  Docking  of  Small  Molecule  Inhibitors 

The  primary  challenge  in  this  report  is  the  computational  discovery  of  small 
molecule  inhibitors.  There  are  tens  of  millions  of  commercially  available  small 
compounds  (molecular  weight:  <500  Daltons)  that  can  be  potentially 
experimentally  tested  to  reduce  HydA  activity  in  bacteria.  Computational 
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prescreening  not  only  greatly  improves  the  cost/benefit  analysis  by  an  order  of  at 
least  1,000,  but  it  also  reduces  the  likelihood  of  false  positives  from  assay  signal 
detection  error  and  unintended  binding  sites. 

We  used  an  automated  computational  pipeline  developed  in-house  that  employs 
Autodock  Vina5  to  dock  small  molecules  to  a  homologous  structure  of  the  C.  acet. 
HydA  from  Clostridium  pasteurianum  (C.  past/,  protein  data  hank  [PDB]  id: 
1FEH1)  (Fig.  2).  Our  group  has  successfully  used  similar  screening  protocols  to 
find  small  molecule  inhibitors  of  host/virus  and  virus/virus  interactions  in  Ebola 
and  Marburg  viruses6'9  and  an  ATP  (adenosine  triphosphate)-powered  secretion 
system  in  plague  bacteria.10  Two  million  commercially  available  compounds  from 
the  ZINC-processed  ChemDiv  catalog11  were  screened  in  2  separate  sites  on  the 
protein  surface.  The  docked  molecules  were  scored  by  consensus  of  the  Vina 
scoring  function  and  the  third-party  DrugScore  DSX  function.  The  top  100 
molecules  from  each  of  the  2  docking  sites  were  visualized  and  10  molecules 
from  each  screen  were  selected  for  purchase.  See  the  Appendix  for  the  structures 
of  the  ZINC  molecules  tested. 


Fig.  2  Atomic  structure  of  closely  related  C.  past,  iron  (Fe)-HydA  protein4 


3.3  Computational  Docking  of  Peptide  Inhibitors 

Docking  of  biological  moieties  such  as  peptide-protein  complexes  has  long  been 
recognized  as  presenting  a  challenge  because  of  a  variety  of  issues  including 

IT 

extreme  flexibility  of  both  components,  the  transferability  of  the  scoring 
function  (the  expression  used  to  assess  the  energy  of  the  system),  and  the  difficult 
necessity  of  a  blind  approach  to  determining  a  complex  where  little  is  known 
beyond  the  structure  of  the  protein  and  the  sequence  of  the  peptide  inhibitor.  Our 
group  has  successfully  explored  these  issues  in  recent  work  developing  peptide 
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binders  to  Anthrax  Protective  Antigen  and  Anthrax  Lethal  Factor,14  in  which  it 
was  demonstrated  that  a  combination  of  molecular  dynamics  and  standard 
docking  protocols  helps  mitigate  these  challenges.15,16 

Inhibition  of  the  C.  acet.  Fe-HydA  has  largely  been  studied  from  the  perspective 
of  small  molecule  inhibitors  such  as  trinitrotoluene  (TNT).  Development  of 
peptide-based  inhibitors  requires  combining  information  from  the  small  molecule 
inhibitors  developed  as  outlined  previously  with  infonnation  obtained  from 
interactions  of  HydA  with  proteins  such  as  ferredoxin  (Fd),  the  physiological 
electron  donor  of  Fe-HydA  in  C.  acet.  In  this  case,  an  additional  level  of 
complexity  is  introduced  because  of  the  absence  of  an  experimentally  derived 
structure  of  either  C.  acet.  HydA  or  Fd.  Both  proteins  are  sufficiently  similar  to 
their  counterparts  in  C.  past,  to  allow  development  of  suitable  homology  models. 
A  model  of  the  docked  HydA:Fd  complex  was  then  developed  and  analyzed.  A 
suitable  peptide  moiety  was  pared  down  from  the  Fd  component  to  provide  a 
reasonable  starting  guess  at  a  peptide  inhibitor  for  the  HydA.  The  peptide  was 
built,  simulated  with  molecular  dynamics,  docked  to  HydA,  and  analyzed  in  the 
context  of  the  successful  small  molecule  inhibitors  determined  in  the  first  part  of 
this  project. 

3.4  Development  and  Testing  of  Gas  Chromatography  Assay 

A  new  protocol  has  been  established  to  test  the  effects  that  potential  HydA 
inhibitors  may  have  on  gas  production  in  C.  acet.  fennentations.  Typically,  gas 
production  is  monitored  at  15 -min  intervals  throughout  the  fennentation  using  gas 
chromatography.  The  minimum  volume  for  running  our  bioreactors  is  normally 
0.5  1.  Thus  it  is  not  practical,  given  limited  availability  (<10  mg)  of  several  of  the 
compounds,  to  screen  using  this  system.  To  circumvent  this  obstacle,  we 
developed  a  protocol  to  perform  these  experiments  in  sealed  serum  vials  and  take 
endpoint  measurements  of  the  headspace  with  a  micro  gas  chromatograph.  For 
this  protocol,  fennentations  were  started  in  an  anaerobic  chamber  containing  an 
atmosphere  of  100%  nitrogen  gas  (No).  An  additional  anaerobic  chamber  was  set 
up  specifically  for  these  experiments,  as  the  atmosphere  of  the  chamber  used  for 
routine  growth  contains  No,  carbon  dioxide  (C02),  and  hydrogen  gas  (H2).  The 
concentrations  of  C02  and  H2  were  measured  after  the  fermentation  was 
complete.  By  comparing  H2:C02  ratios  with  control  fennentations  it  can  be 
detennined  if  a  test  compound  inhibits  the  HydA. 

To  verify  the  method,  we  tested  TNT  since  it  has  been  demonstrated  to  inhibit 
hydrogen  production  from  C.  acet.  fermentations  in  our  laboratory.  Preliminary 
experiments  indicated  that  we  were  able  to  measure  a  30%  reduction  in  hydrogen 
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production  from  fermentations  containing  500-liM  TNT  compared  with  the 
negative  control.  TNT  was  used  as  a  positive  control  in  all  HydA  inhibitor  screens 
to  ensure  that  our  system  was  continuing  to  perform  correctly. 

4.  Results  and  Discussion 


4.1  Structural  Model  of  the  Hydrogenase 

We  tried  to  generate  an  accurate  model  of  HydA  by  perfonning  a  careful  analysis 
on  a  wide  array  of  homology  models  (Fig.  3).  Multiple  comparative  modeling 
algorithms  were  tested  for  the  construction  of  a  C.  acet.  HydA  homology  model, 
including  SWISS-MODEL,18  CPHmodels,19  and  I-TASSER.20  Agreement  among 
these  models  was  excellent,  as  well  as  agreement  with  the  C.  past,  crystal 
structure  (PDB  id:  1FEH1).  Alignment  and  comparison  with  the  1FEH  structure 
gave  a  backbone  root  mean  squared  deviation  (RMSD)  of  3.44  A  for  the 
CPHmodels  algorithm,  3.30  A  for  I-TASSER,  and  3.32  A  for  SWISS-MODEL.  In 
addition  to  an  excellent  RMSD,  I-TASSER  was  the  only  one  of  these  algorithms 
that  did  not  truncate  the  C.  acet.  HydA  structure  to  agree  with  the  somewhat 
shorter  C.  past,  protein.  I-TASSER  was  also  used  to  develop  a  homology  model 
for  the  C.  acet.  Fd  based  upon  the  C.  past,  crystal  structure  (PDB  ID:  1CLF). 


Fig.  3  Overlay  of  3  homology  models  of  the  C.  acet.  HydA  protein  generated  by  1) 
CPHModels  (gray),  2)  SWISS-MODEL  (yellow),  and  3)  I-TASSER  (orange)  with  the 
experimental  structure  (PDB  ID:  1FEH,  blue) 

Unfortunately,  we  found  that  protein-protein  docking  using  comparative  models 
of  this  resolution  can  lead  to  anomalous  results  whereby  side  chains  are  not  in  the 
right  place  or,  worse,  entangled  with  each  other.  For  the  computational  screening 
of  compounds,  we  chose  to  use  the  C.  past,  structures  of  Fd  and  HydA  as  is.  The 
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risk  was  that  we  could  have  found  the  right  molecules  for  the  wrong  protein. 
However,  we  observed  good  sequence  similarity  between  C.  past,  and  C.  acet. 
HydA  in  the  docking  target  regions  we  have  chosen. 

4.2  Protein-Protein  Docking  to  Determine  Sites  for  Inhibitors 

To  determine  the  sites  to  computationally  dock  to,  we  used  2  well-regarded 
protein-protein  docking  programs,  ZDOCK  and  Patch  Dock,'  to  model  how  the 
2  proteins  connected  to  each  other  ’  ’  (Fig.  4).  With  a  significant  amount  of 
molecular  surface  area  in  the  HydA,  we  found  at  least  3  putative  docking 
locations  for  the  much  smaller  Fd  protein.  Of  these  3  locations,  2  provided  a  close 
approach  between  the  2  Fe-sulfur  (S)  clusters  of  Fd  (blue)  and  the  Fe-S  cluster 
“wire”  in  the  HydA  (red).  Close  approach  is  required  for  electron  transfer  from 
the  Fd  “shuttle”  to  the  HydA  “port”.  The  electrons  absorbed  by  the  HydA  wire  are 
used  to  drive  the  chemical  reaction,  2H+  +  Fd  (reduced)  — >  2H2  +  Fd  (oxidized), 
which  results  in  more  H2  output  and  less  overall  butanol  production  in  the 
bacteria. 


Fig.  4  Putative  docking  site  1  of  Fd  (blue)  to  Fe-only  HydA  (red);  Fe-S  clusters  are  green 
and  yellow 

4.3  Computational  Small  Molecule  Docking  of  the  2  Sites 


For  each  of  the  2  possible  docking  sites,  we  computed  a  midpoint  coordinate  for  a 
cubic  box  with  a  side  length  of  20  A.  We  then  extracted  the  database  of 
3 -dimensional  structures  of  compounds  from  the  commercial  vendor,  ChemDiv 
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(San  Diego,  CA),  from  the  ZINC  database.11  With  roughly  1.2  million  entries,  we 
filtered  out  molecules  with  predicted  log  partition  coefficients  >4.0  to  ensure 
that  we  did  not  consider  “greasy”  molecules  that  tend  to  stick  to  many  proteins 
and  could  cause  unintended  inhibition  of  other  bacterial  processes.  The  remaining 
approximately  1  million  molecules  were  computationally  docked  with  Autodock 
Vina5  within  the  specified  boxes.  The  molecules  were  free  to  flex  around  their 
torsional  angles,  and  often  several  docking  structures  were  generated  for  each 
compound.  Twenty  structures  for  each  of  the  approximately  1  million  candidate 
molecule  were  scored  with  the  third-party  DrugScore  program.  The  scores  from 
Vina  and  Drugscore  were  nonnalized  and  summed  to  yield  a  Z-score  that  was 
used  to  rank  the  best  molecules.  We  visualized  the  top  50  scoring  molecules  for 
each  docking  site,  and  chose  the  structures  that  looked  best  by  expert  opinion.  The 
top  10  compounds  from  each  site  were  labeled  by  their  original  ChemDiv  catalog 
number  and  the  lists  were  sent  to  the  SEDD  experimental  team. 

4.4  Computational  Peptide  Docking  of  the  Entire  Protein 

Analysis  of  the  Patchdock-generated  HydA:Fd  complex  showed  a  tightly  folded 
single  loop  of  the  Fd  consisting  of  residues  21-33  (sequence:  NAISQGDSIFVID) 
to  be  situated  within  the  site  1  crease  of  the  HydA  in  close  proximity  to  the  Fe2S2 
cluster.  This  was  taken  as  an  excellent  starting  point,  and  a  peptide  of  that 
sequence  was  built,  solvated,  minimized,  and  simulated  at  nonnal  pressure  and 
temperature  of  1  atm  and  300  K  for  40  ns  with  a  time  step  of  2  fs  using  the 
NAMD  software."  It  was  apparent  during  the  course  of  the  simulation  that  this 
single-fold  configuration  of  the  peptide  is  moderately  stable,  as  the  strand 
unfolded  and  refolded  multiple  times.  Multiple  configurations  reflecting  both 
folded  and  unfolded  states  were  selected  from  the  simulation  trajectory  and 
docked  to  HydA  using  Patchdock.  This  represents  a  global  dock  on  the  entire 
HydA  structure. 

Ranked  analysis  of  the  results  showed  a  “hotspot”  of  folded  configurations  within 
the  site  1  crease,  as  expected,  while  an  additional  hotspot  for  unfolded 
configurations  was  seen  along  the  side  leading  to  site  2  (Fig.  5).  Electrostatic 
analysis  with  the  APBS  software  of  Baker  et  al.  highlights  the  co-localization  of 
electropositive  (blue)  residues  in  each  of  these  sites,  allowing  favorable 
interaction  with  peptide  electronegative  residues  at  the  end  and  center.  The  folded 
configuration  of  the  peptide  binding  at  site  1  allows  the  segregation  of  a  cluster  of 
hydrophobic  residues  at  one  end  of  the  peptide  onto  the  protein  interface  while 
allowing  solvent  exposure  of  hydrophilic  residues  clustered  on  the  upper  fold  of 
the  peptide.  The  unfolded  peptide  configurations  fit  more  neatly  into  the  extended 
patch  of  electropositive  surface  residues  adjacent  to  site  2.  There  are  interesting 
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structural  and  chemical  similarities  between  the  folded  peptide  configuration  and 
small  molecule  inhibitors  determined  in  the  first  part  of  the  project.  Both  sets  of 
inhibitors  are  approximately  the  same  size  and  shape  for  fitting  the  site  1  patch 
shown  in  Fig.  5,  occupying  a  slightly  bent  but  overall  rectangular  volume  of 
roughly  20  A  on  the  long  axis.  The  chemical  character  of,  for  example,  small 
molecule  candidate  inhibitor  ZINC21938598  is  analogous  to  the  lower  fold  of  the 
peptide,  consisting  of  electronegative  groups  at  each  end  suitable  for  binding  to 
the  protein  with  a  string  of  hydrophobic  aromatic  rings  in  the  middle.  We  propose 
that  the  folded  configuration  of  the  peptide  is  optimum  for  binding,  as  it 
aggregates  hydrophobic  residues  on  one  side  and  hydrophilic  on  the  other.  The 
logical  next  extensions  of  this  work  are  to  explore  mutations  to  further  stabilize 
this  folded  configuration  and  to  make  single-point  mutations  on  the  binding  side 
of  the  fold  to  further  enhance  binding. 


Fig.  5  Peptide  docking  on  HydA  shows  hotspots  at  site  1  (central  right)  and  side  entrance  to 
site  2  (lower  left).  Electrostatic  surface  of  protein  was  generated  with  APBS  software.27 
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4.5  Testing  for  Hydrogenase  Inhibition 


In  the  first  round  (Fig.  6),  we  tested  all  20  compounds  as  previously  described  at 
500  uM  and  found  that  compounds  H,  I,  J,  and  T  appeared  to  reduce  hydrogen 
production  in  these  fermentations.  We  also  noted  that  compounds  F  and  N  killed 
the  bacteria,  which  led  us  to  perfonn  a  second  round  of  testing  at  lower 
concentrations  (Fig.  7)  to  ensure  the  reason  for  lethality  was  not  due  to  too  strong 
of  a  reduction  of  HydA  activity. 
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Fig.  6  H2/C02  ratio  after  incubation  with  small  molecules 


9 


■  H2  (%) 

■  CO 2  (%) 

■  H2/C02 


1.600 
1.400 
1.200 
1.000 
0.800 
0.600 
0.400 
0.200 
0.000 

y  y  y  >  > 

.  k  .  i'  —i'  _  1/  x  1/  x  ^  ^  ✓  F  ^  v  F  </V  A  A  A 

^  «  «  O  O  ^  ^  V  V  s  s 

c#  c#  °  ° 


1 

cn 

. 

■ 

1 

n 

m 

1 

yyyyyyyyyyyyyy  y  y  y  y  y  y 
o°Ao>J?i’  'W  <<N  <<N  <A  <A  <A  <<N  V  ,. V  ,. V  ,. V  V  -V  ■  v 


^  <£>  <5X'  <3' 
(J'  o2. 


<<<<<<  <c  r  V 


..y\.y 

,  V  .  V 


<< y  s  ’ 


■  H2  (%) 

■  C02  (%) 

■  H2/C02 


Fig.  7  Representative  results  from  tests  for  HydA  inhibition,  (top)  Gas  concentrations  in 
culture  headspace  after  6  h  of  incubation  with  test  compounds.  DMSO  was  used  a  solvent  for 
the  test  compounds  and  serves  as  a  negative  control.  TNT  is  a  positive  control  for  the  assay, 
(bottom)  Gas  concentrations  form  cultures  receiving  dilutions  of  compound  F. 


After  a  second  round  of  testing  (Fig.  7;  data  not  shown),  we  concluded  that  none 
of  the  molecules  selected  by  computational  screening  provided  a  statistically 
significant  reduction  in  H2  production.  In  the  future,  we  will  look  at  a  few  of  the 
following  possibilities: 

•  Crystallization  and  X-ray  structure  determination  of  the  actual  C.  acet. 
Fe-HydA  protein.  This  protein  is  of  sufficient  interest  to  the  hydrogen  fuel 
production  community  to  be  worthwhile  to  characterize  more  carefully. 

•  Experiments  to  determine  the  exact  binding  site  of  Fd  to  the  HydA.  Not 
only  would  this  provide  a  site  for  computational  docking,  but  it  would  also 
help  researchers  understand  how  electrons  flow  from  the  Fd  to  the  HydA. 

•  Computational  screening  of  the  new  protein  structure  or  computational 
screening  of  a  sufficiently  accurate  homology  model. 

•  Selection  and  purchase  of  peptide  fragments  resembling  the  Fd  protein  as 
potential  peptide  inhibitors  of  the  HydA. 
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5.  Conclusions 


A  new  experimental  assay  for  determining  H2  production  in  ultra-small  bacterial 
cultures  was  developed  and  verified  with  TNT.  In  addition,  small  molecule 
inhibitors  of  C.  acet.  Fe-HydA  were  computationally  selected.  Moreover,  a 
peptide  inhibitor  has  been  designed  computationally  from  the  Fd:HydA  complex 
that  shares  attributes  of  the  selected  small  molecules.  Several  of  the  putative  small 
molecule  inhibitors  were  found  to  reduce  H2  production  moderately,  but  not 
significantly.  Other  molecules  actually  reduced  cell  viability.  Future  directions 
were  suggested. 


11 


6.  References 


1.  Servinsky  MD,  Germane  KL,  Liu  S,  Kiel  JT,  Clark  AM,  Shankar  J,  Sund  CJ. 
Arabinose  is  metabolized  via  a  phosphoketolase  pathway  in  Clostridium 
acetobutylicum  ATCC  824.  J  Ind  Microbiol  Biotechnol.  2012;39:1859-1867. 

2.  Cai  X,  Servinsky  M,  Kiel  J,  Sund  C,  Bennett  GN.  Analysis  of  redox 
responses  during  TNT  transformation  by  Clostridium  acetobutylicum  ATCC 
824  and  mutants  exhibiting  altered  metabolism.  Appl  Microbiol  Biotechnol. 
20 13;97(10):456 1-4563. 

3.  Lutke-Eversloh  T,  Bahl  H.  Metabolic  engineering  of  Clostridium 
acetobutylicum:  recent  advances  to  improve  butanol  production.  Curr  Opin 
Biotechnol.  2011;22:634-647. 

4.  Peters  JW,  Lanzilotta  WN,  Lemon  BJ,  Seefeldt  LC.  X-ray  crystal  structure  of 
the  Fe-only  hydrogenase  (CpI)  from  Clostridium  pasteurianum  to 

1.8  angstrom  resolution.  Science.  1998;282:1853-1858. 

5.  Trott  O,  Olson  AJ.  AutoDock  Vina:  improving  the  speed  and  accuracy  of 
docking  with  a  new  scoring  function,  efficient  optimization,  and 
multithreading.  J  Comput  Chem.  2010;31:455-461. 

6.  Liu  Y,  Lee  MS,  Olson  MA,  Harty  RN.  Bimolecular  complementation  to 
visualize  lilo virus  VP40-host  complexes  in  live  mammalian  cells:  toward  the 
identification  of  budding  inhibitors.  Adv  Virol.  201 1  ;20 1 1:341816. 

7.  Han  Z,  Lu  J,  Liu  Y,  Davis  B,  Lee  MS,  Olson  MA,  Ruthel  G,  Freedman  BD, 
Schnell  MJ,  Wrobel  JE,  Reitz  AB,  Harty  RN.  Small-molecule  probes 
targeting  the  viral  PPxY-host  Nedd4  interface  block  egress  of  a  broad  range 
of  RNA  viruses.  J  Virol.  2014;88:7294-7306. 

8.  Lu  J,  Han  Z,  Liu  Y,  Liu  W,  Lee  MS,  Olson  MA,  Ruthel  G,  Freedman  BD, 
Harty  RN.  A  host-oriented  inhibitor  of  Junin  Argentine  hemorrhagic  fever 
virus  egress.  J  Virol.  2014;88:4736-4743. 

9.  Brown  CS,  Lee  MS,  Leung  DW,  Wang  T,  Xu  W,  Luthra  P,  Anantpadma  M, 
Shabman  RS,  Melito  LM,  MacMillan  KS,  Borek  DM,  Otwinowski  Z, 
Ramanan  P,  Stubbs  AJ,  Peterson  DS,  Binning  JM,  Tonelli  M,  Olson  MA, 
Davey  RA,  Ready  JM,  Basler  CF,  Amarasinghe  GK.  In  silico  derived  small 
molecules  bind  the  filovirus  VP35  protein  and  inhibit  its  polymerase  cofactor 
activity.  J  Mol  Biol.  2014;426:2045-2058. 


12 


10.  Swietnicki  W,  Carmany  D,  Retford  M,  Guelta  M,  Dorsey  R,  Bozue  J,  Lee 
MS,  Olson  MA.  Identification  of  small-molecule  inhibitors  of  Yersinia  pestis 
type  III  secretion  system  YscN  ATPase.  PLoS  One.  2011  ;6(5):e  1 9716. 

1 1 .  Irwin  JJ,  Shoichet  BK.  ZINC  -  a  free  database  of  commercially  available 
compounds  for  virtual  screening.  J  Chem  Inf  Model.  2005;45:177-182. 

12.  Velec  HF,  Gohlke  H,  Klebe  G.  DrugScore(CSD)-knowledge -based  scoring 
function  derived  from  small  molecule  crystal  data  with  superior  recognition 
rate  of  near-native  ligand  poses  and  better  affinity  prediction.  J  Med  Chem. 
2005;48:6296-6303. 

13.  London  N,  Raveh  B,  Schueler-Furman  O.  Peptide  docking  and  structure- 
based  characterization  of  peptide  binding:  from  knowledge  to  know-how. 
Curr  Opin  Struct  Biol.  2013;23:894-902. 

14.  Sellers  M,  Hurley  M.  In  silico  design  of  smart  binders  to  anthrax  PA. 
CBRNE  Sensing  XIII.  In:  Fountain  AW,  editor;  Proceedings  of  Chemical, 
Biological,  Radiological,  Nuclear,  and  Explosives  (CBRNE)  Sensing  XIII; 
2012  Apr  24-27;  Baltimore,  MD.  Bellingham  (WA):  International  Society 
for  Optics  and  Photonics;  p.  835807-835809. 

15.  Sellers  MS,  Hurley  MM.  XPairlt  docking  protocol  for  peptide  docking  and 
analysis  of  molecular  simulations.  Submitted  to  Molecular  Simulations; 
2015. 

16.  Sellers  MS,  Hurley  MM.  US  Anny  Research  Laboratory  (ARL)  XPairlt 
Simulator  for  Peptide  Docking  and  Analysis.  Aberdeen  Proving  Ground 
(MD):  Army  Research  Laboratory  (US);  2014  Jul.  Report  No.:  ARL-TR- 
6999.  Also  available  at:  http://www.arl.armv.mil/arlreports/2014/ARL-TR 
-6999.pdf. 

17.  Demuez  M,  Coumac  L,  Guerrini  O,  Soucaille  P,  Girbal  L.  Complete  activity 
profile  of  Clostridium  acetobutylicum  [FeFe]-HydA  and  kinetic  parameters 
for  endogenous  redox  partners.  FEMS  Microbiol  Lett.  2007;275: 1 13-121. 

18.  Biasini  M,  Bienert  S,  Waterhouse  A,  Arnold  K,  Studer  G,  Schmidt  T, 

Kiefer  F,  Cassarino  TG,  Bertoni  M,  Bordoli  L,  Schwede  T.  SWISS-MODEL: 
modeling  protein  tertiary  and  quaternary  structure  using  evolutionary 
information.  Nucleic  Acids  Res.  2014;42:W252-258. 

19.  Nielsen  M,  Lundegaard  C,  Lund  O,  Petersen,  TN.  CPHmodels-3.0  -  remote 
homology  modeling  using  structure-guided  sequence  profiles.  Nucleic  Acids 
Res.  2010;38:W576-581. 


13 


20.  Zhang,  Y.  I-TASSER  server  for  protein  3D  structure  prediction.  BMC 
Bioinformatics.  2008;9:40. 

21.  Chen  R,  Li  L,  Weng  Z.  ZDOCK:  an  initial-stage  protein-docking  algorithm. 
Proteins.  2003;52:80-87. 

22.  Schneidman-Duhovny  D,  Inbar  Y,  Nussinov  R,  Wolfson  HJ.  PatchDock  and 
symmdock:  servers  for  rigid  and  symmetric  docking.  Nucleic  Acids  Res. 
2005;33:W363-367. 

23.  Chang  CH,  King  PW,  Ghirardi  ML,  Kim  K.  Atomic  resolution  modeling  of 
the  Fd:[FeFe]  HydA  complex  from  Chlamydomonas  reinhardtii.  Biophys  J. 
2007;93:3034-3045. 

24.  McCullagh  M,  Voth  GA.  Unraveling  the  role  of  the  protein  environment  for 
[FeFe]-HydA:  a  new  application  of  coarse-graining.  J  Phys  Chem  B. 
2013;117:4062-4071. 

25.  Cheng  T,  Zhao  Y,  Li  X,  Lin  F,  Xu  Y,  Zhang  X,  Li  Y,  Wang  R,  Lai  L. 
Computation  of  octanol-water  partition  coefficients  by  guiding  an  additive 
model  with  knowledge.  J  Chem  Inf  Model.  2007;47:2140-2148. 

26.  Phillips  JC,  Braun  R,  Wang  W,  Gumbart  J,  Tajkhorshid  E,  Villa  E,  Chipot  C, 
Skeel  RD,  Kale  L,  Schulten  K.  Scalable  molecular  dynamics  with  NAMD. 

J  Comput  Chem.  2005;26:1781-1802. 

27.  Baker  NA,  Sept  D,  Joseph  S,  Holst  MJ,  McCammon  JA.  Electrostatics  of 
nanosystems:  application  to  microtubules  and  the  ribosome.  Proc  Natl  Acad 
Sci  USA.  2001;98:10037-10041. 


14 


7.  Transitions 


At  this  time,  no  publication  or  presentations  of  this  work  beyond  this  report  are 
planned.  However,  we  fully  expect  scientific  results  and  technological 
improvements  from  this  work  will  inform  future  endeavors  in  our  respective 
laboratories. 
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Appendix.  Small  Molecules  Experimentally  Tested 


This  appendix  appears  in  its  original  form,  without  editorial  change. 
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Table  A-l  List  of  small  molecules  experimentally  tested  in  this  work 


ID  ZINC  Database  ID 


Structure 
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Table  A-l  List  of  small  molecules  experimentally  tested  in  this  work  (continued) 


ID  ZINC  ID 


Structure 


F  ZINC21938598 


I  ZINC09429018 


19 


Table  A-l  List  of  small  molecules  experimentally  tested  in  this  work  (continued) 


ID  ZINC  ID 


Structure 
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Table  A-l  List  of  small  molecules  experimentally  tested  in  this  work  (continued) 


ID  ZINC  ID 


Structure 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


ARL 

C02 

C.  acet. 

CEBP 

C.  past. 

Fd 

Fe 

HydA 

H2 

n2 

PDB 

RMSD 

S 

SEDD 

TNT 


FIS  Army  Research  Faboratory 
carbon  dioxide 
Clostridia  acetobutylicum 
commonly  expressed  bacterial  protein 
Clostridia  pasteurianum 
ferredoxin  (iron-sulfur  protein) 
iron 

Fe-Fe  Hydro  genase  A 
hydrogen  gas 
nitrogen  gas 

Protein  Data  Bank  (online  repository  of  protein  structures) 

root  mean  squared  deviation 

sulfur 

ARE  Sensors  and  Electron  Devices  Directorate 
trinitrotoluene 
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